Derivative expansions of real-time thermal 

effective actions 

Maria Asprouli, Victor Galan-Gonzalez 

Imperial College 
Theoretical Physics Group 
Blackett Laboratory 
London SW7 2BZ 

February 1, 2008 

Abstract 

In this work we use a generalised real-time path formalism with 
properly regularised propagators based on Le Bellac and Mabilat Q 
and calculate the effective potential and the higher order derivative 
terms of the effective action in the case of real scalar fields at finite 
temperature. We consider time-dependent fields in thermal equilib- 
rium and concentrate on the quadratic part of the expanded effective 
action which has been associated with problems of non-analyticity at 
the zero limits of the four external momenta at finite temperature. We 
derive the effective potential and we explicitly show its independence 
of the initial time of the system when we include both paths of our 
time contour. We also derive the second derivative in the field term 
and recover the Real Time (RTF) and the Imaginary Time Formalism 
(ITF) and show that the divergences associated with the former are 
cancelled as long as we set the regulators zero in the end. Using an 
alternative method we write the field in its Taylor series form and we 
first derive RTF and ITF in the appropriate limits, check the analytic- 
ity properties in each case and do the actual time derivative expansion 
of the field up to second order in the end. We agree with our previous 
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results and discuss an interesting term which arises in this expansion. 
Finally we discuss the initial time-dependence of the quadratic part 
of the effective action before the expansion of the field as well as of 
the individual terms after the expansion. 

1 Introduction 

The interest in the amalgamation of field theory and statistical mechanics 
arose from the realisation that many problems encountered experimentally 
and theoretically in particle physics have many-body aspects. For this rea- 
son, zero-temperature quantum field theory was reformulated by generalising 
the usual time-ordered products of operators to the ordered products along 
a path in a complex time-plane [[|. The choice of the path gives rise to 
different formalisms but all theories should give the same physical answers. 
Although the various path-ordered finite temperature field theory formalisms 
such as Real Time Formalism (RTF) including the closed-time approach || 
and Imaginary Time Formalism (ITF) [|J should give the same physics, there 
has been serious discussion about their exact equivalence. 

In this paper we will tackle a problem in RTF which consists in the 
occurrence of pathologies associated with singularities, arising in diagrams 
with self-energy insertions or in some effective potential calculations. This 
problem appears when products of delta functions with the same argument 
are involved and creates non-analyticities in the effective action at finite 
temperature, thus making it ill-defined. 

The main interest for developing an effective formalism which describes 
the finite temperature field theory comes from the need to tackle important 
problems in phase transitions, which have played a crucial role in the early 
evolution of the universe. The significance and observable quantities of a 
specific transition depend on its detailed nature and its order. A reason for a 
well-defined effective action comes from the fact that it represents the quan- 
tum corrections which in general might be of extreme importance in defining 
or changing the order of a transition. For example, analytic analysis || sug- 
gests that the electroweak phase transition is first order because of quantum 
corrections from gauge bosons while non-perturbative lattice simulations of 
high temperature electroweak theory suggests that this is only true if the 
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Higgs is lighter than 70Gev f|. The inclusion of higher derivative terms in 
the derivative expansion of the effective action in a first order transition is of 
great importance in cases such as the derivation of the rate of the sphaleron 
fluctuations. These configurations have been used to explain the observed 
baryon asymmetry of the universe M. 

Moreover, although the effective potential can give the approximate crit- 
ical temperature of a given transition, it is not adequate for answering ques- 
tions concerning the departure of the field from equilibrium occurring during 
dynamical cooling near and below the critical temperature T c . The effective 
potential describes static properties and therefore it is not an appropriate 
tool for studying the dynamical behaviour of a wide class of field theory 
models considered in inflationary scenarios. 

The standard method for estimating the quantum corrections is to first 
integrate out quantum fluctuations about a constant background. This gives 
an effective potential for cf) which is then used in the equations of motion 
determining (p(x,r) 0. Integrating out fluctuations about a general inho- 
mogeneous configuration gives the full effective action which includes the 
higher derivative terms. In the language of quantum field theory at finite 
temperature the effective potential r is given by r = j3F, where F is the 
minimum of the free energy at which the system lies in the case of local 
thermal equilibrium. If the ensemble averages of the matter fields are homo- 
geneous and static, then the free energy is given by the finite temperature 
effective potential ||. 

In analogy to quantum mechanics, the decay rate of an unstable config- 
uration (fif with energy £f is given by |L0| 

T = ~lm[£f) = "2Im[ lim ^\n(<p f \e' HT / h \cf) f )} 

where H is the Hamiltonian and the matrix element can be described as a 
functional integral 

(<l>f\e- BT " i \<j> f ) = N jvct>e- s ^l h 

in Euclidean time. Here <fi is subject to the condition <p(T/2) = <ft(—T/2) = 
<f)f and 5* denotes the Euclidean action. Evaluating the functional integral 
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to one loop, we expand S(<f>) about a solution of the equations of motion, 0, 
and keeping only terms quadratic in the fluctuations £0 = — 0, we obtain 

NjV(j)e- s ^ /n ~ iVe- s( * /7i) [det(-«9^ + ^"(0))]^ 
= exp[-3 eff {^/ti)] 

where S e ff is the effective action. If we expand S e ff about a constant 0, 
i.e. in powers of momentum about a point with zero external momenta, in 
position space and zero temperature this reads 

S eff {<j>) = J d 4 x[ -K//(0) + \z{(j>)d^4> + 0((d^) 4 ) ] 

where we have made use of T = Lorentz properties. For constant only 
the effective potential term survives. 

Although such an expansion up to the second derivative has been per- 



formed at zero temperature for scalar and Dirac field theories |fLl]j , there are 
difficulties arising in the equivalent expansion at finite temperature. Das and 
Hott [|T^] find a non-analyticity in the two-point functions involving the tem- 



perature dependent term of the quadratic part of the effective action. In this 
spirit, if the derivative expansion breaks down at finite temperature, the def- 
inition of the effective potential might not be unique. This non-analyticity, 
in the case of the vacuum polarisation for a scalar field coupled to a classical 
external field, manifests itself in the difference between the order of the lim- 
its (p — > 0, p — > 0) and (p — > 0,p — > 0) of the external momenta, the first 
relating to the electric screening mass of the photon and the second to the 



plasma frequency of the particular field theory under consideration [13|. This 



non-commutativity appears in hot QCD |]14 |, self interacting scalars [15], [16 



and gauge theories with chiral fermions HI7| . In ITF, setting p M = first and 
performing the mode sum gives the same result as taking the limit po — > 
first and the limit p —>■ afterwards. In RTF extra Feynman rules have 
been imposed to explain this difference in the two limits [0. The problem 
is neither due to subtleties in the use of Feynman parametrisation at finite 



temperature Jl9[ , nor to the infinite number of possible extensions of po to the 



imaginary axis and its analytic continuation to the complex plane ||16||. The 
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lack of analyticity and the infrared divergences occurring in the definition 
of the effective action at finite temperature show the need for an effective 
field theory formalism from which RTF and ITF rules can be derived easily 
In the next chapter we will describe a method dealing with these problems. 
We will calculate the two-propagator contribution (bubble diagram) to the 
second derivative term in the effective action using Le Bellac and Mabilat's 
generalised real-time path formulation with properly regularised propagators 



2 The method 

In Le Bellac and Mabilat's approach [|IJ they derive Feynman rules that 
take explicitly into account the vertical part of the contour and recover the 
RTF in the case of diagrams with at least one finite external line and the 
ITF in the case of vacuum fluctuations. They keep the regulators of the 
propagators when they find problematic products of delta functions with the 
same argument and show that they can use RTF when no such problems 
arise. They claim that the contribution of the vertical part of the contour 
lies in the cancellation of the t{ dependent terms of the horizontal part since 
the whole result should be U and £/ independent due to the KMS condition 
of the propagators. We will now describe in detail this method, which we 
will use throughout this paper. 



2.1 Outline of the method 



The specific approach uses the Mills |Kj mixed representation of the prop- 
agators for a free scalar field, with t defined in the generalised time path C, 
starting at t{ and ending at tj — if] of Fig.l. The propagator is written as 

D c (t,k) = J ^ e -*»*[e c (t) +n(ko)]p(k ,-k) (1) 

where 6 c (t) is a contour 9 function, n(ko) is the Bose-Einstein distribution 
function given by 
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Figure 1: The integration contour C in the complex t-plane. 

and p(ko, k) is the (temperature independent) two-point spectral function 
given by 

p(k ,k) = 27re(k )5(k 2 -ujl) (3) 
where e(ko) is the sign function and 



OA 



m 



(4) 



However, the propagators need to be regularised because eventually we 
want to take the Fourier transform in time by taking the limits ti — > — oo, tf — > 
+oo (to get energy conservation) and this is ill-defined since the integrands 
are linear combinations of complex exponentials. For this reason we write 
the S distribution in its regularised form 



S(k =F v k ) 



1 



1 



2in k =f — ie k =f oj^ + ie 
and thus p(k , k) in Eq. (§) can be written as 

i rs 



p(fc ,k) 



E 



2u k r f^ ±l k - su k + ier 



The regularised propagator can be written as 

D c R (t, k) = D>(t, k)6 c (t) + D<(t, k)0 c (-t) 



(5) 



(6) 



(7) 
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and obeys the KMS condition 



D>(t - t^k) = D<(t,k) (8) 

Momentum integration in the complex ko-pl&ne will give for D^(t,k), 
with t defined in the region —f3 < Imt < 

D>{t) = ^-Y.\ 0(e) + n(u k - ie*) ]e^^ - £ X v e~^ (9) 
where s = sign[i?e(t) ], £ is the regulator and X,, is the sum 

x = y rs = %& r ,ujk ^ 

where u v = ^jp = 2irrjT denotes the Matsubara frequencies. This term arises 
from the residues of the distribution function when we integrate ko in the 
complex plane. We will discuss its contribution later. 

2.2 The bubble term 

The rest of the paper will be the calculation of the bubble diagram, which 
is nothing else but the product of two propagators. In order to justify its 
relevance, we will briefly mention where it comes from. We consider the two 
scalar field theory described by the Lagrangian 

£[<M = \d^V - l -m 2 r, 2 - X -g^ 2 + £ (11) 

where Cq denotes the free Lagrangian for <fi. If we integrate out the 77-field 
fluctuations and use a one-loop approximation we find that the generating 
functional can be expressed as 

Z = ( V(j)e lSom+iS 'cff m (12) 
J c 

where S' e ff is given by 



S' eff [<P] = ilHn[l - gA c (x,x')<f>(x')} (13) 
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and A c (x,x f ) is the propagator for the r\ field. Expanding the logarithm we 

get 

oo 

S'eff[<l>] = Y, S eff (14) 
p=l 

where p denotes the number of the propagators. In this expression we will 
concentrate on the first non-1 
the expansion and is given by 



(2) 

concentrate on the first non-local term which is the quadratic part S e Jf of 



(2) -ig 2 Tr(A c (x,x')(f)(x')A c (x',x)(/)(x)) 
S eff = ^ (15) 



One can permute the order of the elements inside the trace using the identity 

<f>(x)A c (k) = A c {k + p)<p{x) (16) 



which is equivalent to the Taylor expansion of Fraser 21] for moving mo- 
mentum operators to the left of functions depending on x, when we identify 
p^ = —idp. Thus the quadratic part of S' e fA(j)] can be rewritten as 



7 (2) 
'e// 



--J d A x J d*x'(j)(x)iB(p,l3)(j)(x') (17) 
with iB(p,ft) being the bubble term given in terms of the propagators as 



iB(p, (3) = g 2 [ ^ / A c (k, m)A c (k + p, m) (18) 

(2) 

Separating the time dependence which interests us at finite temperature, S e jj 
is written 

,2 fti _i P d s 



(2) %g- r^p f <ry 



r (i 3 k fti-i-i 3 

/ -—3 / dt'A c (k; t, t')A c {k + p; t', t)0(p, (19) 
J [Ztc) Jti 



and the fields are periodic over the time path ^ — i(3\. This derivative 
expansion of the bubble term is well established at zero temperature. At finite 
temperature this is not so, since we will be expanding our theory around 
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an ill-defined point. This can be seen for example if we look at the An 
component of the propagator of our theory which is given by 



A(k,m) = — 2i7rn(k )5(k 2 - m 2 ) (20) 

k z — m l + is 

Substituting the propagator in Eq. (|17D, the temperature-dependent real part 
of the quadratic thermal effective action is given by the following expression 
which is nonanalytic at the zero four-momentum limit (T2j 



with 



2 



oc 



Re(S eff ^[ ( p]) = -^j d 4 x^ dkf{k)cj>{x)\cj>{x) (21) 



_ kn(u)Re(\nR) 
tW- o;(-V 2 )V2 

where uj 2 = k 2 + m 2 and 

R= H (d 2 - V 2 + 2isud + 2irk{-V 2 ) 1/2 ) r (23) 

r,s=+l,-l 

This result suggests that the derivative expansion breaks down at finite tem- 
perature, due to the problematic product of the two delta functions contained 
in the bubble. If the derivative expansion of the effective action is not rigor- 
ously possible, then the definition of an effective potential, the lowest order 
term in such an expansion, is not unique. This would have consequences 
in any kind of study concerning symmetry breaking and restoration, unless 
there is a formal way to overcome these pathologies and have a well defined 
derivative expansion. 



2.3 Effective potential term 

We use Le Bellac and Mabilat's formulation [[]J to prove explicitly the tj- 
independence of the effective potential in the case of two propagators and 
one external time to (first term of the bubble diagram) when adding both 
contributions from the horizontal and vertical path. Since we are interested 
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in time dependent fields, we will concentrate on the time integral of the 
bubble term in Eq. (|19|). This is given by 

= / D c R (t - tjD&to - t x )dt x (24) 

where C will be in the region [U,to] and [to,ti] for the horizontal path and 
[ti,ti — i/3] for the vertical one as shown in Fig.l. Using the representation 
for the propagators of Eq. (H), the integrations over the two paths give G H 
for the horizontal and Gy for the vertical one |L| 

G R = _ J±«<L-1 * )[ i _ e — (ei+e 2 )(t -u) 2au-t ) i (25) 

{2uuf y iu(e 1 + e 2 ) + 2£ ;[ 

Gy = ; 6 — L -)l e -Me i+ e 2 )(t -U) 2S(U-t o)] (26) 

v {2uf y tu(s 1 + e 2 ) + 2^ n 

where 

Any = E [ ^( ±£ +n(u- i& x ) ] [ 0(±e 2 ) +n(u- i£e 2 ) ] (27) 



are ^-independent coefficients. The KMS relation which reads as 

[ 0(e) + n{uj - i£e) } = e^e^l 9{-e) + n(uj - i£e) } (28) 

gives for the coefficients Ah and Ay 

A H = e^+^e'^Av (29) 

Adding the tj-dependent terms of both paths Eq. ( ^5f) and Eq. ( p6|) and using 
the KMS condition of Eq. (p8|), we find that they cancel. Therefore, it is 
essential that we add the vertical contribution to ensure the ij-mdependence 
of our result. The remaining part of the sum is given by 

G R = G H + Gy = - J-p J£ H + (30) 

and using the definitions of Ah and Ay from Eq. (B7I) and identities of the 6 1 
functions, the sum is written as 



rR+rR- 1 j (1 + 2n(^))i^ + 2i£ 2 ri(u) \ i/3n(u)(l + n(u)) 
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In Eq. (pnp the first term is the result of Eq. fl30| ) for E\ = 82 and the second 
term is the result for E\ + 82 = 0. Taking the limit of the regulator £ to zero 
at the end, we have 

G g + G « = ^(l + 2„M) + ^MiL±£M (32) 

The previous result agrees completely with the ITF result for the effective 
potential [p^l , which proves the consistency of our theory to this order. Now 



we examine the cases of taking different limits for the £j and the regulators 
and try to explain their physical meaning. 

1. We take the limit ti — * — 00 (keeping the regulators finite) which should 
give us the real time formalism. Since the total sum is tj-independent, 
the limit of ti — > — 00 is already given by Eq. in which the regu- 
lator can be taken to zero since there is no need for it any more after 
the limit has been performed. We notice from Eq. ( |26| ) that keeping 
the regulators finite the vertical part vanishes in this limit, recovering 
thus the RTF and the total sum is being given by the tj-independent 
contribution of the horizontal part. In the case of unregularised propa- 
gators, the vertical contribution contains a tj-independent term of the 
form 

£_L n (o,)(l + n (uj))2Tt5(k 2 - m 2 ) 
i 2cu 

which in the regularised approach is hidden in the two horizontal parts 
of the contour, as seen in the last part of Eq. ([31~1). We see that, to this 
order, the regularised formalism is dealing with the pathologies of the 
problematic delta functions, recovering RTF in the appropriate limit. 

2. Now keeping ti finite, we take the zero limits of the regulators in dif- 
ferent orders and find 

+ s 2 = 0, i - 0) = G R H (i - 0, s 1 + s 2 = 0) = 

and the only contribution comes from the vertical part in the limit 
81 + 82 = 0, £ — > which is the second term of Eq. (j3ll) 

Inn c= ^Hd + "("» 

£->0,ei+£ 2 =0 2u 2 
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This limit is one part of the full effective potential of Eq. ( p2[ ) as ex- 
pected since it only corresponds to the limit of equal and opposite 
energies ui for the two propagators (ei + £2 = 0). 

2.4 The second derivative term 

Now the same formalism will be used for the derivation of the second deriva- 
tive term of the effective action in our 1-loop case, where only time- dependent 
fields are considered. Because now the sum will also contain terms polyno- 
mial in (ti — to) as well as exponential ones, the equivalence with the RTF 
and ITF is less straightforward. We will show that the RTF limit can be 
extracted in this case without problems of divergences as long as we keep the 
regulators finite and take them to zero after the limit has been done. The 
second derivative term will look like 

r^ 2) = f dhD c R {t , ti)(* - hYD^h, t ) (33) 

where we have omitted the ^d 2 (f>(to) factor of this term in the expansion. 
Using the definition of Eq. (fy), T 2 can be written as 

^ = Ic dtl I I Vcfati) +n(k ))(9 c (t 1 ,t ) +n(h)) 

x(ti -t ) 2 e~ i(fco ~ fcl ~ i£)(i °^ 1 V(^o)p(^i) (34) 

where the regulator e is used so that the limit of ti — > —00 can be taken 
without problems. In the end it will be set to zero. The other two regulators 
E\ and So in the delta functions of p(k )p(ki) make sure that no problems 
appear in the equal energy (mass) case k = k\ = w. 

Now we first perform the dt\ integration and then "absorb" the (ti — to) 2 
term by differentiating the result with respect to ko- If we name the time 
integrals over the two paths as Ic, we then write 

/ dhih - to) 2 e e {t lf ; ) e -(*o-*i--)(*o-*i) = (35) 
Jc oko 

Substituting this formula into our general expression Eq. fl34|), the contribu- 
tions from the different paths can now be written 

rg) = J ^ (n(ko) _ „(fc 1 ))p(*b)p(* 1 )i a ^/H (36) 
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r gj = J ^^ [n{ko){n{kl) + l)]^)^ 2 J^I y (37) 

In the 9 2 term vanishes due to the opposite sign of its time arguments. 
The n 2 term vanishes due to the cancellation between the two horizontal 

(2) 

paths. In T v one of the ^-functions always vanishes due to the choice of to 
on the horizontal path. 

The time integrations Ih, Iv over the two paths give 

(\ — e — i(ko-ki— ie)(to—ti)\ 

Ih = i : ; ; 

^0 — M — ^ 

(~\ — pP{ko-k\—ie)\ 
j _ i e -i(ko-ki-i£)(t -ti) x ^ 1 



ko — ki — ie 

The analytical calculation of the different path contributions is quite com- 
plicated since it involves first and second order residues and therefore deriva- 
tives of the distribution functions. We performed the momentum integrations 
and then took the same limits of our variable At = t~ t and of the regulators 
as before to check the consistency of our method for the second derivative 
term. 

1. We took the At — > — oo limit keeping the regulators finite. In the total 
sum the At-dependence appears in terms like At n e eAt and At n (n = 
0, 1, 2). These terms could cause divergences in the At — > — oo limit but 
they disappear once we include the vertical part in our calculation. Our 
result is independent of the order in which the regulators are taken to 
zero in the end and is given by a finite term coming from the horizontal 
part 

lim r f2)_ » (l + 2n(o;)) 

This term looks like the first order term of the effective potential divided 
by u 2 , as it can be seen from Eq. (|32]), which is sensible since it is 
essentially the first correction due to the second derivative. 

2. Our second limit is tj = to,e = in order to try to recover the ITF 
result (tj = to = finite and e is not needed any more since U is finite). 
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This proved to be also independent of the order of the zero limits of 
the regulators. We obtained 

r< 2 >(ti = to,e = 0;ei = 0;eo = 0) = 



i 



(2n(u) + 1) - (3u{2n{u){n{u) + 1) + 1) 



+ ^V(2, i(l ,) + l) + 4gV " ( " )( " ( " ) + 1) ] 

3 

which is consistent with the derivation of the second derivative term in 
the ITF formalism [E31. 



3 Alternative method 

Another possible way of performing our calculation is to consider the full 
Taylor series of the field but do the actual expansion and study the individual 
terms in the end. We will generalise our method considering different energies 
(uj and VL) in the delta functions of Eq. @ for each propagator of the bubble 
term. In this way we will check the analyticity limits of the full derivative 
term by taking the limits Q — > ±u ( V — > 0) and <9 t — > in different orders 
at the end of the calculation. 

The expanded field can be written as 



<f>(ti 



co -r Qn 



n=0 



<f>(t) 



t=t 



t=t 



The full derivative term of the field inserted between the two propagators, 
which is the last time integral in Eq. looks like 



JC *=*o 



(38) 



This term acts as an energy-shift by —idt in the exponentials of the propa- 
gators making the time-integrals over the paths Ic of Eq. (|35|) look like 



c 



dt 1 9 c (t 1 ,t Q )e 



-i(fco— fci 



-i(e+9 t ))(to-ti) _ p 
— 1 C 



(39) 
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Performing the time integration for the horizontal and vertical path as before, 
we get 

_ g-i(fco-fci-i(e+3t))(to-t»)) 

I>H = l h - h - i(e + d t ) 

and 

/i _ P(k -ki-i(£+d t ))\ 
j/ _ ^ e -i(fco-fci-i(e+9t))(*o-*i) x v 1 c >_ 

V k Q -k 1 -i{e + d t ) 

The energy integration gives us the full bubble term as a sum over the two 
paths written in terms of At = U — t 

T (B) = T (B) + r (B) 

with 

-p(B) _ ^ ^- Ko )n(»- Kl ) r ( e ^+n-»( e o+ e i))_i)(i_ e -^A f) i 
ff ~ ^ 4jfi X[ A J [ } 



and 

(B) _ in(w - ie )n(tl - ie^ e^e-^^e^-l) 



where 

A = cj + f2 — + e — e — d t ) 

Now we can check the analyticity of our result keeping At finite. We expand 
the distribution functions and take the limits of our regulators to zero (we 
can do that since we keep At finite). If we take the limits Q — > ±uj and 
dt — > 0, we get finite and independent of the order of the limits results. The 
full derivative expansion of the bubble term, therefore, is analytical in this 
limit and is 

r<f> (n - ±., a, - o) = , (2 » M + (42 ) 

for the horizontal case and 

rS? ) ( n - to, a, - o) = ^M + iW*^) + + D (43) 



2a; 3 u> 2 

for the vertical one 
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1. Now we consider the limit U = to in Eq. ( f4~2"|) and Eq. ((4"3|). This gives 

r? ] = (44) 

and 

This is exactly the result for the effective potential using the ITF for- 
malism, as expected since it is the zeroth time and space-derivative 
term, when U = t . It also agrees with our previous result in Eq. (|32]) 
of the effective potential after we set the regulators to zero. (In Eq. 
(|32|) the result differs by a factor of 1/2 due to the fact that we have 
initially considered same energies u for the propagators and this corre- 
sponds to half of the result of the effective potential of Eq. (fi5D when 
different energies are assumed). 

2. Now we take the limit At — > — oo and we will check the analyticity 
again. In this case only the first At-independent part of T H survives 
and taking the regulators to zero after the limit has been performed, 
we get 

r<*>(At - -oo) = £ (46 ) 



The analyticity check for Eq. ([46]) gives finite but different results for 
different orders of performing the limits (Q — > ±u, d t — > 0). We found 
that performing the time-derivative (dt — > 0) limit first and the spatial- 
derivative (Q — » ±u) afterwards, we had the usual effective potential 
term of Eq. ([45|) 

lim r < B > = f ( 2 ^H + 1 ) + . Pn(u)(l + n(u)) 

but reversing the order of the limits gave us only the first term of our 
previous result 

Um rW = i (2»M + D 

n->±u>,dt-+Q 2u 3 
We see that although we don't have divergence problems in taking the 
limits in both orders, approaching the zero from the space-derivative 
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first seems to produce only part of the full result in agreement with the 
result of Evans using ITF [p2]] . Now we take only the spatial derivative 



to zero in the At — > — oo case of Eq. Q4q) which gives 

Umr <B ) = 2i (^M+0 

If we now expand our result in powers of the time-derivative d t , we get 
the zeroth order term of our analyticity check and a second order term 
of the form 

i(2nH + l) 2 



This is exactly the second order time-derivative term derived in this 
limit using our previous method in section 2.4. We have to note that 
in the case of At — > — oo, there is no term linear in dt. 



Now we perform the same expansion in powers of the time-derivative 
up to the second order but for a general finite At, for both horizontal and 
vertical paths in Eq. ([40"D and Eq. (41) and take the limits (Q — > ±u) to get 



the effective potential and the higher derivative terms. We identify the terms 
as follows 

1. d® term 

r ^ = -i (2n( ^ 3 +1) (e 2 ^ + e ~ 2wAt - 2) (47) 



u 2 



8} term 



F l H = (M^) + 1) [(e 2^At _ e -2^A t) _ 2icu(At)(e 2i " At + e ^ At )} (49) 
8oj 
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1 



F ^ = -7-7[(2nM + l)(e 



leu 4 



2iuiAt _ e ~2iu>At\ 



- 2/3u((n(u) + l)V^ A * - n 2 (uj)e- 2lu)At ) - A(3 2 u 2 n(uj)(n(uj) + 1)] 
+ ^ [(2n(w) + 1) (e 2lujAt + e~ 2iuM ) 

+ A(3un(uj)(n(uj) + l)](At) (50) 

We notice the existence of a non-zero <9t-dependent term unlike the zero 
temperature case where such a term vanishes. This could be related to 
the loss of Lorentz invariance in the finite temperature case and could 
be interpreted as an energy shift. The existence of such a linear term 
might be of great physical importance in the study of time-dependent 
systems. Such a term did not exist in the expansion for the At infinite 
case, where only zero and second order terms in the time-derivative sur- 
vived. This makes sense since in the infinite time limit any interaction 
with the heat bath which gives rise to such linear terms will have been 
damped. Mathematically this term could arise due to the shape of the 
time contour, which in the finite At case is non-symmetric. However 
this is not the case for the zero-temperature situation or the non zero 
temperature one in the infinite At limit where the symmetry of the 
contour will make any time integration of odd terms in the derivative 
expansion to vanish. In the At = case this term is equal to 

r i _ -»/?r° 

2 

where T° is the effective potential term given by Eq. QT7D and Eq. (fi8|) 
in the At = limit. 

3. d\ term 

F ^ = Z(2 ^ 5 + ^ K e2iujAt + e ~ 2iujAt - 2 ) - 2tu(At)(e 2iujAt - e ~ 2iujAt ) 

- 2iu 2 (At) 2 (e 2iujAt + e- 2iujAt )] (51) 



T 2 v = r 2 yi + r 2 y2 (52) 
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with 

T 2 V1 = --^[(2n( W ) + l)(e 2,wAl + e~ 2iujAt ) - 2{3uo((n(uj) + l) 2 e 2iujAt 

+ n 2 (u)e- 2 ^ At ) 

+ 2/3V((n(w) + lf e 2iwAt - n 2 (uj)e- 2iujAt ) 

+ ^Vn( U )(n( W ) + l)] (53) 

and 

T 2 V2 = —A(2n(uj) + l)(e 2lujAt - e - 2iujAt ) - 2/3w((n(w) + lf e 2iu}At 

- n 2 {u)e~ 2iuAt ) - A(3 2 u 2 n{uj){n{uj) + l)](At) 



+ -^-[(2 n ( W ) + l)(e 2 ^ + e 
+ A(3ujn(uo) (n(w) + 1)]( At) 2 (54) 

If we take the ij = to limit of our second derivative term, we recover 
our previous derivation of the same term in section 2.4. 

In our calculation we have omitted the contribution of the X v term of Eq. 
10). This term which arises from the residue of the distribution function 
vanishes since it is proportional to the regulator. In the finite At case the 
regulators are set to zero before any limit is taken while in the infinite At 
case they are set to zero once the infinity limit has been performed. In both 
cases this term does not contribute. 



2iujAt , -2iu>At^ 



4 The initial time dependence 

In this section we will treat the initial time-dependence of our problem in a 
rather more formal way. 

Based on Le Bellac and Mabilat's proof of the tj-fndependence of a regu- 
larised Green function [0], we will prove the same for our effective potential 
term. Our tj-dependent integrals in this case are 

f° dhGnih, t )G R (t u t ) + ^ dhGniti, t Q )G R {t u t ) (55) 
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Differentiating the first term with respect to U we get 

—Gr(U, to)Gft(ti, t ) 
Repeating for the second term we now get 

G>(t i -tf3,t )G>{t l -t(3,t Q ) 
Using the KMS condition for thermal equilibrium 

G>(t - i/3) = G<(t) 

we see that these terms cancel. If we repeat the same method for the higher 
derivative terms of the field of the bubble case, we get a non-zero result which 
shows explicitly the tj-dependence of these terms. The same analysis for the 
second derivative term gives 

G<(U, t ) [-P 2 - 2if3(U - t )]G<(U, to) (56) 

If we generalise in the case of the m-th derivative term, the ti dependence of 
the derivative with respect to t^ will have the form 

G$(Mo)E II ^r^H^fe - to) m - k ]G< R {UM) (57) 
k=l j=o K - 

We see that the individual terms of the expanded field are clearly t ^-dependent 
even in the case of tj = t , where the highest order /5-term survives in the 
previous sum. This is somehow expected since a truncated expansion of the 
field makes it no longer periodic. On the other hand if we have a periodic 
field (p(t) in equilibrium, the derivative with respect to ti discussed earlier 
will give 

- G<(t u toMtJGZiU, to) + G>(U - *o)0(*i - ifiG^U - i(3, t ) (58) 

which is zero for periodic fields and regularised propagators obeying the KMS 
condition. 
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5 Conclusions and possible applications 



We found that using this closed time path formalism we can avoid the 
pathologies in the RTF and derive the ITF limit as well, both in the case of 
the effective potential and in the second derivative correction of the bubble 
diagram. The fact that we can cancel the divergence in our effective action 
using our prescription, hence providing us with a formalism which allows us 
to compute quantum corrections to the effective potential is the key point 
of our paper. However we found that the inclusion of the vertical path and 
the careful treatment of the regulators are essential for the cancellations to 
happen. 

We showed a general way to compute higher derivative terms in the bub- 
ble and derived the complete bubble term. We checked its analyticity for 
finite and infinite time differences At and found different limits in the second 
case. The non zero, linearly dependent on the time derivative of the field 
term found in the finite At case, may be related to the loss of translational 
invariance at finite temperature. The physical meaning of such a term and 
in particular its sign and whether it is complex or real may be important in 
the study of time-dependent systems. Gribosky and Holstein |l3j do not find 
such a linear term in their expansion of derivatives of the field. Motivated 
by the use of Feynman parametrisation at zero temperature \1M they calcu- 
late the vacuum polarisation diagram using ITF but extending to continuous 
Po first and evaluating the mode sum afterwards. They compare their re- 
sult with Dittrich's |23] background field method of calculating the effective 
Langrangian at finite temperature and in both cases there is no linear term 
unlike our case which appears for any finite At. 

The extension of our calculation to higher derivative terms and to space- 
dependent fields will give a full effective action whose importance in field 
theory was discussed earlier. Our calculation may be performed for higher 
order diagrams in the expansion of the one-loop effective action, but this is 
beyond the scope of this paper. We have considered a two real scalar field 
theory, but we could in principle use our method in different models, such as 
a Yukawa or a gauge theory or even consider systems with time-dependent 
parameters. The possibility of evaluating quantum corrections can be more 
directly applied to phase transitions, where they may indicate us something 
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about the order of the transition. 
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